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Natural convection in air-filled rectangular cavities inclined with respect to gravity, so that the heated wall is fac- 


ing upwards, is studied numerically under the assumption of two-dimensional laminar flow. A computational code 


based on the SIMPLE-C algorithm is used for the solution of the system of the mass, momentum and energy 


transfer governing equations. Simulations are performed for height-to-width aspect ratios of the enclosure from 


0.25 to 8, Rayleigh numbers based on the length of the heated and cooled walls from 10° to 10’, and tilting angles 


of the enclosure from 0° to 75°. The existence of an optimal tilting angle is confirmed for any investigated confi- 
p g ang y g 


guration, at a location that increases as the Rayleigh number is decreased, and the height-to-width aspect ratio of 


the cavity are increased, unless the value of the Rayleigh number is that corresponding to the onset of convection 


or just higher. Dimensionless correlating equations are developed to predict the optimal tilting angle and the heat 


transfer performance of the enclosure. 


Keywords: Natural convection in air, Tilted cavity, Differential heating, Numerical analysis, Dimensionless cor- 


relations 


Introduction 


Natural convection inside differentially-heated rec- 
tangular enclosures filled with air has been widely stu- 
died in the past decades, due to its relevance to many 
engineering applications, such as heat transfer in build- 
ings, solar energy collection, and heat removal in mi- 
cro-electronics, to name a few. 

The main body of the research conducted on this topic 
deals with situations wherein the enclosure is either ver- 
tical or horizontal, which means that buoyancy is induced 
by imposing a heating either from the side, in examining 
conventional convection, or from below, in studying 
thermal instabilities. However, a number of studies re- 
lated to the effects of the cavity inclination with respect 
to gravity on the enhancement or deterioration of the heat 
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transfer performance of the enclosure are also available 
in the literature, demonstrating the existence of an op- 
timal tilting angle of the enclosure for maximum heat 
transfer [1-19]. A summary of these publications, subdi- 
vided according to the research method, 1.e., experimen- 
tal or numerical, is presented in Tables 1 and 2, respec- 
tively, in which indications on the investigated ranges of 
the tilting angle of the enclosure with respect to the grav- 
ity vector, ọ (where ọ = 0° corresponds to the vertical 
configuration in which the cavity is differentially heated 
at sides, and positive tilting angles denote configurations 
with the heated wall facing upwards), the height-to-width 
aspect ratio of the enclosure, A, and the Rayleigh number 
based on the cavity height, Ray, are enumerated. Indica- 
tions on the optimal tilting angle of the enclosure, which 
our interest is focused on, are also reported. Indeed, 
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Nomenclature 
A height-to-width aspect ratio of the enclosure 
g gravity vector, m/s? 
g gravitational acceleration, m/s” 
H height of the enclosure, m 
k thermal conductivity, W/(m K) 
Nu Nusselt number 
P dimensionless pressure 
Pr Prandtl number 
Q heat transfer rate, W 
q heat flux, W/m? 
Ra Rayleigh number based on the cavity height 
T dimensionless temperature 
t temperature, K 
U,V  x,y-wise dimensionless velocity component 
V dimensionless velocity vector 
W width of the enclosure, m 


X, Y dimensionless Cartesian coordinates 
xX, y Cartesian coordinates, m 
Greek symbols 
a thermal diffusivity, m?/s 
B coefficient of volumetric thermal expansion, 1/K 
v kinematic viscosity, m?/s 
p mass density, kg/m’ 
T dimensionless time 
tilting angle of the enclosure with respect to 
? gravity, deg 
y dimensionless stream function 
Subscripts 
av average 


cooled wall, at the temperature of the cooled wall 
heated wall, at the temperature of the heated wall 


opt optimal value 


Table 1 Summary of experimental studies performed on natural convection in air-filled inclined cavities 


Year Author(s) Ray H/L 0) opt 

1975 Ozoe et al. [1] 6.6 x 10°—2.2 x 10’ 8.4, 15.5 -90° — 90° ~ 10° 

1976 Hollands et al. [2] 2x 10°- 1.1 x 10° 48 30° = 75° 

1982 Elshirbiny [3] 1.3 x 105- 2.4 x 10" 5-110 0° — 90° 

1983 Elshirbiny [4] 1.0 x 10°— 1.5 x 10" 10, 80 0° — 30° 

1985 Schinkel and Hoogendorn [5] 2 x 10°—2.7 x 108 1-11 0° -70° 15° — 30° 

1988 Karayiannis and Tarasuk [6] 6.0 x 10°— 1.8 x 10° 6.7 — 16.7 30° — 75° 

1989 Hamady et al. [5] 1 x10- 1x106 1 0° — 180° 10° — 20° 

2007 Bairi et al. [8] 4.2 x 10° -4.2 x 10’ 0.75 0° — 360° 60° — 90° 

2008 Bairi [9] 1x 10'-1 x10" 1 0° — 360° 60° — 90° 

2012 Corvaro et al. [10] 1.2 x 10°-—4.6 x 10° 1 0° — 90° 0° 
Table 2 Summary of numerical studies performed on natural convection in air-filled inclined cavities 

Year Author(s) Ray H/L 0) opt 

1983 Elshirbiny [4] 1.0 x 10°—1.5 x 10!° 10, 80 0%- 30° 

1993 Kuyper et al. [11] 1x 10-1 x 10” 1 0° — 180° 

1996 Soong et al. [12] 5.4 x 10*— 1.3 x 10° 1,3,4 0%- 90° 0° — 30° 

2001 Delgado-Buscalioni et al. [13] 1x 10°-8 x 10° 0.01 — 0.25 0° — 90° 0° — 20° 

2006 Baez and Nicolas[14] 1 x 10°- 5.12 x 10° 1,8 0°— 330° 

2010 Nor Azwadi et al. [15] 1x 10° 1 20° 90° 10° — 30° 

2011 Munir et al. [16] 5x 10° 1 0°— 180° ~ 10° 

2012 Vivek et al. [17] 1 x 10f— 1 x 10° 0.2-1 0°— 180° ~ 15° 

2014 Chang [18] 1 x 10°-5.12 x 108 1-8 0°— 180° 0° — 20° 

2016 Williamson et al. [19] 1x 10*—1 x 108 1 0°— 80° 15°— 25° 


although some results obtained by the cited authors are 
contradictory, in line of principle the optimal tilting angle 
was found to increase as both the Rayleigh number and 
the aspect ratio of the enclosure were decreased. Howev- 
er, none of the cited authors developed a correlation for 


the prediction of the optimal tilting angle and the corres- 
ponding heat transfer rate across the cavity. Moreover, 
the physics behind the existence of the optimal tilting 
angle and its location in relation with the value of the 
independent variables used as controlling parameters was 
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generally left out of the discussion of the results. 

Framed in this general background, a comprehensive 
study of natural convection in air-filled, differentially- 
heated rectangular cavities inclined with respect to grav- 
ity, so that the heated wall is facing upwards, is per- 
formed numerically under the assumption of two-dimen- 
sional laminar flow, for ranges of the height-to-width 
aspect ratio of the enclosure from 0.25 to 8, the Rayleigh 
number based on the length of the heated and cooled 
walls from 10° to 10’, and the tilting angle of the enclo- 
sure with respect to gravity from 0° to 75°, with the main 
scope to investigate the basic heat and momentum trans- 
fer features and develop accurate correlations for pre- 
dicting the optimal tilting angle and the heat transfer 
performance of the enclosure, which seems useful for 
thermal engineering design purposes. 


Mathematical formulation 


An air-filled rectangular enclosure of height H and 
width W is considered. The enclosure is differentially 
heated at two opposite walls, which are kept at uniform 
temperatures t, and t,, while the other pair of walls are 
assumed to be perfectly insulated. The enclosure is in- 
clined at an angle with respect to the gravity vector, 
where = 0° denotes the vertical configuration in which 
the cavity is differentially heated at sides. A sketch of the 
enclosure is displayed in Fig. 1, where the reference Car- 
tesian coordinate system (x, y), and the thermal states of 
the boundary walls, are also represented. The resulting 
buoyancy-induced fluid flow is considered to be 
two-dimensional, laminar and incompressible, with con- 
stant fluid properties and negligible viscous dissipation 
and pressure work. The buoyancy effects on momentum 
transfer are taken into account through the customary 
Boussinesq approximation. 


Fig.1 Sketch of the geometry and coordinate system 


Once these assumptions are employed in the conserva- 
tion equations of mass, momentum and energy, the fol- 
lowing set of governing dimensionless equations is ob- 
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tained: 
V-V=0 (1) 
N 4 (v-vyve-ve+v27v- 2478 (2) 
OT Pr g 
oT VVT -Lyr (3) 
OT r 


where V is the dimensionless velocity vector having x- 
wise and y-wise dimensionless velocity components U 
and V normalized by v/H, t is the dimensionless time 
normalized by H’/v, P is the dimensionless pressure 
normalized by pv’/H’, T is the dimensionless temperature 
excess over the uniform temperature of the cooled wall 
normalized by the temperature difference between the 
heated and cooled walls (th — te), g is the gravity vector, 
Pr = v/a is the Prandtl number, and Ra is the Rayleigh 
number defined as 


pa = Solin =f) (4) 
av 
where B, v and a are the coefficient of volumetric 
thermal expansion, the kinematic viscosity and the ther- 
mal diffusivity of the fluid, respectively. 
Another parameter which enters into this study is the 
height-to-width aspect ratio of the enclosure 
H 
ae (5) 
The assigned boundary conditions are: (a) T= 1 and V 
= 0 at the heated wall; (b) T = 0 and V = 0 at the cooled 
wall; and (c) ôT/ðn = 0 and V = 0 at the adiabatic walls, 
where n denotes the normal to the wall. 
The initial conditions assumed throughout the enclo- 
sure are fluid at rest, i. e., V = 0, and uniform tempera- 
ture T= 0. 


Computational procedure 


The system of the governing equations defined by Eqs. 
(1)-(3), combined with the boundary and initial condi- 
tions stated earlier, is solved through a control-volume 
formulation of the finite-difference method. The pres- 
sure-velocity coupling is handled using the SIMPLE-C 
algorithm introduced by Van Doormaal and Raithby [20], 
which is a more implicit variant of the SIMPLE algo- 
rithm developed by Patankar and Spalding [21]. Convec- 
tive terms are approximated by the QUICK discretization 
scheme proposed by Leonard [22]. A second-order 
backward scheme is employed for time integration. Both 
space discretization and time stepping are chosen uni- 
form. Starting from the assigned initial fields of the de- 
pendent variables across the cavity, at each time-step the 
system of discretized algebraic governing equations is 
solved iteratively by way of a line-by-line application of 
the Thomas algorithm. A standard under-relaxation tech- 
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nique is enforced in all steps of the computational pro- 
cedure to ensure adequate convergence. Within each time- 
step, the spatial numerical solution of the velocity and tem- 
perature fields is considered to be converged when the 
maximum absolute values of the mass source, as well as 
the relative changes of the dependent variables at any grid- 
node between two consecutive iterations, are smaller than 
the pre-specified values of 10° and 10°, respectively. As 
time passes, the dynamic behavior of the system ana- 
lyzed is followed by plotting the phase trajectories of the 
dependent variables T, U, and V, at a number of assigned 
spatial locations, i.e., the distributions of the local time- 
derivatives of the dependent variables versus the varia- 
bles themselves with time as a parameter, and continuo- 
usly monitoring the time distributions of the heat transfer 
rates at the heated and cooled walls, as well as their rela- 
tive difference, so as to visualize the path of the system 
toward the steady-state. Time-integration is stopped when 
the relative changes of the time-derivatives of the depen- 
dent variables at any grid node and the relative difference 
between the heat transfer rates added to the fluid by the 
heated wall and withdrawn from the fluid by the cooled 
wall between two consecutive time steps are smaller than 
the pre-set values of 10” and 10°, respectively. 

Once steady-state is reached, the local Nusselt num- 
bers at the heated and cooled walls, Nu, and Nu,, are 
calculated by the following expressions: 


Niis qH _ oT (6) 
k@,—t,)  OX|y-9 

Nu, - Se a i 
k(te=t)  OX|y asa 


where qr and q, are the heat fluxes at the heated and 
cooled walls, k is the thermal conductivity of the fluid, 
and X and Y are the dimensionless Cartesian coordinates 
normalized by H. Both temperature gradients in Eqs. (6) 
and (7) are evaluated by a second-order temperature pro- 
file embracing the wall-node and the two adjacent 
fluid-nodes. The corresponding average Nusselt numbers 


(Nup)ay and (Nu,),, are then calculated as: 
1 


Q, ôT 
Np jpa — i=j ay 8 
(Nuna = E J al; (8) 
1 
Q ôT 
Njen eel ay (9) 
i k(te—tp) la X=1/A 


0 

where Qh and Q. are the heat transfer rates added to 
the fluid by the heated sidewall and withdrawn from the 
fluid by the cooled sidewall. Both integrals in Eqs. (8) 
and (9) are computed numerically by means of the trape- 
zoidal rule. Of course, since at steady-state the heat 
transfer rates at both the thermally active walls are the 
same, i.e., Qn = —Q, = Q, the average Nusselt numbers 
coincide with what can be assumed to be the average 
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Nusselt number of the enclosure Nu, = (Nup)ay = (Nug)ay- 
Numerical tests related to the dependence of the re- 
sults obtained on the mesh spacing and time stepping 
have been performed for several combinations of the 
three controlling parameters, namely A, Ra, and ọ. The 
mesh-spacing and time-step used for computations are 
chosen in such a way that further refinements do not 
produce noticeable modifications either in the heat trans- 
fer rates or in the flow fields. In particular, the percentage 
change of Nu, as well as that of the maximum absolute 
value of the dimensionless stream function throughout 
the enclosure |W|max, Where y is defined as usual by U = 
Ow / OY and V = —Ow / OX, must be smaller than the pre- 
established accuracy value of 1%. The typical number of 
nodal points and the time-step used for computations lie 
in the ranges between 60 x 60 and 120 x 480, and between 
10° and 10%, respectively. Selected results of the grid 
and time-step sensitivity analysis are listed in Tables 3-5. 
Finally, with the scope to validate the numerical code 
used for the present study, three different tests have been 
executed. In the first test, the steady-state solutions ob- 
tained for an air-filled, vertical square cavity differen- 
tially heated at sides with adiabatic top and bottom walls 
have been compared with the benchmark solutions of de 
Vahl Davis [23], Mahdi and Kinney [24], Hortman et al. 
[25], and Wan et al. [26]. In the second test, the average 
Nusselt numbers computed numerically at steady-state 
for vertical rectangular cavities using a Prandtl number 
Pr = 0.7 and Rayleigh numbers Ra = 10°—5x10’, have 
been compared with the usually recommended Ber- 
kovsky—Polevikov correlation, see, e.g., Bejan [27] and 
Incropera et al. [28]. In the third test, the average Nusselt 
numbers computed numerically for inclined rectangular 
cavities have been compared with the results of Bairi [9]. 
A good degree of agreement between our numerical solu- 
tions and the literature data was achieved in any valida- 
tion test carried out, as e.g. shown for the third test in Fig. 
2. More details on the code validation can be found in a 
previous paper authored by Corcione and Habib [29]. 
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Fig. 2 Comparison between present numerical results and the 
results of Bairi [9] for a square enclosure 
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Table 3 Time-step sensitivity analysis for square cavity 
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Ra A [deg] Mesh size At [s] Nun % TE liraz % 
10° 1 0° 120 x 120 10° 9.47 - 25.53 - 
10° 9.11 -3.85 24.56 -3.80 
10* 8.93 -1.96 24.09 -1.94 
10° 8.93 0.01 24.06 -0.11 
Table 4 Grid sensitivity analysis for inclined square cavity 
Art (s) Ra A ọ[deg] Mesh size Nun % [P| max % 
10° 10° 1 0° 60 x 60 4.60 - 13.95 - 
80 x 80 4.55 -1.22 13.76 -1.38 
100 x 100 4.54 -0.12 13.71 -0.37 
30° 60 x 60 4.72 - 21.97 - 
80 x 80 4.64 -1.58 21.59 -1.75 
100 x 100 4.64 -0.17 21.54 -0.22 
60° 60 x 60 4.51 - 34.19 - 
80 x 80 4.45 -1.49 33.65 -1.59 
100 x 100 4.43 -0.30 33.49 -0.47 
Table 5 Grid sensitivity analysis for inclined rectangular cavity 
Ar (s) Ra o[deg] A Mesh size Nu, % [P| max % 
10° 10° 0° 1 60 x 60 9.33 - 25.28 - 
80 x 80 9.16 -1.81 24.79 -1.93 
100 x 100 9.13 -0.32 24.68 -0.43 
120 x 120 9.11 -0.23 24.56 -0.49 
0° 0.25 60 x 240 61.92 - 30.62 - 
80 x 320 61.21 -1.15 30.02 -1.98 
100 x 400 61.04 -0.27 29.87 -0.48 
120 x 480 60.88 -0.26 29.72 -0.51 
30° 0.25 80 x 320 304.12 - 101.63 - 
100 x 400 298.22 -1.94 99.50 -2.10 
120 x 480 296.80 -0.48 98.75 -0.75 
140 x 560 296.12 -0.23 98.43 -0.32 


Results and discussion 


Numerical simulations are performed for different 
values of (a) the height-to-width aspect ratio of the cavity, 
A, in the range between 0.25 and 8, (b) the Rayleigh 
number based on the length of the thermally active walls, 
Ra, in the range between 10° and 10’, and (c) the tilting 
angle of the enclosure with respect to the gravity vector, 
~, in the range between 0° and 75°. 


Typical local results are reported in Fig. 3, in which 
equispaced isotherm and streamline contour plots at dif- 
ferent tilting angles are depicted for A= 1 and Ra = 10°. 
In the streamline plots, the contours correspond to equis- 
paced values of the normalized dimensionless stream 
function lwl/ | Wl max in the range between Oand 1.In the 
isotherm plots, the contours correspond to equispaced 
values of the dimensionless temperature in the range be- 
tween 0 and 1. 
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It may be seen that at ọ = 0° the flow field, as ex- 
pected, consists of a single-roll cell that derives from the 
rising of the hot fluid adjacent to the heated wall and its 
descent along the opposite cooled wall, which leads to 
the distinctive temperature distribution featured by a flu- 
id stratification in the core of the cavity. As the tilting 
angle ọ is increased, the vertical path available for the 
acceleration of the heated fluid increases, which results in 
a growth of the motion intensity and a simultaneous 
promotion of the thermal stratification breakdown, with a 
consequent increase of the heat transfer rate. On the other 
hand, as ọ is further increased, the jet of cold fluid that 
moves downwards along the adiabatic wall tends to sep- 
arate from it, with the consequent formation of a stagna- 
tion region on the heated wall downstream the point of 
impingement of the descending fluid, and, therefore, a 
more or less significant decrease of the amount of heat 
transferred across the cavity. This is displayed in Fig. 4, 
where the distributions of the local Nusselt number along 
the heated wall, Nu,, are reported for three different tilt- 
ing angles, demonstrating that the thermal performance at 
= 15° is higher than at ọ = 0° and 60°. Equivalent but 


= 


= 


[Wlmax = 13.71 


Wla = 17.14 
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reversed distributions could be depicted for the opposite 
cooled wall. Notice that the single-roll flow structure is a 
common pattern for all the examined configurations. 
Flow pattern evolutions not too dissimilar from that de- 
scribed above are encountered for other values of A and 
Ra, as e.g. shown in Figs. 5 and 6 for A= 0.25 and Ra = 
10°, and for A= 4 and Ra = 10°, respectively. 

As far as the overall results are concerned, the distri- 
butions of the average Nusselt number Nu plotted versus 
Ra at ọ=0° are displayed in Fig. 7 using the aspect ratio 
of the cavity as a parameter. It is apparent that, as long as 
conduction is the dominant heat transfer mode and then 
the Nusselt number is independent of the Rayleigh num- 
ber, the thermal performance of the cavity increases with 
increasing the aspect ratio, due to the increase of the uni- 
form temperature gradient existing between the heated 
and cooled walls. On the other hand, once the onset of 
motion takes place, which occurs at a critical Rayleigh 
number whose value increases with increasing the aspect 
ratio, the Nusselt number increases with increasing the 
Rayleigh number, tending to become almost independent 
of the aspect ratio when the buoyancy force is high 


—— 


Wimax = 27.83 Wimax = 33.49 


Fig. 3 Isotherm and streamline contours plots for Ra = 10°, A= 1 and ọ = 0° — 60° 
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Fig. 4 Local Nusselt number for Ra = 10°, A= 1 and ọ = 0° — 
60° 
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Fig.5 Isotherm and streamline contours plots for Ra = 10°, A 
= 0.25 and ọ = 0° — 45° 
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enough to give rise to the formation of two boundary 
layers adjacent to the thermally active walls. Actually, 
owing to the smaller amount of heat transferred between 
the fluid streams moving along the adiabatic walls, the 
thermal performance is slightly higher for the tall cavities 
and lower for the shallow cavities. An alternative repre- 
sentation of the results reported in Fig. 7 can be obtained 
by plotting the distributions of the average Nusselt num- 
ber Nu versus A using the Rayleigh number as a parame- 
ter, as shown in Fig. 8, in which the mentioned slight 
increase of the thermal performance occurring with the 
increase of the cavity aspect ratio is pointed out. Trends 
of this same type of those illustrated in both Figs. 7 and 8 
can be observed when the cavity is inclined with respect 
to the gravity vector. 

Thus, in order to clearly highlight in what measure the 
inclination with respect to gravity affects the thermal 
performance of the cavity, the distributions of the ratio 
(Nu/Nup) between the average Nusselt numbers for the 
tilting angle ọ and for the vertical position at same Ray- 
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leigh number and aspect ratio, plotted versus @ using Ra 
as a parameter, are depicted in Figs. 9-11 for A= 0.25, 1, 
and 4, respectively. It is apparent that, owing to the pro- 
gressively larger inclination of the cavity, and the men- 
tioned increase of the motion intensity, the ratio (Nu/Nuo) 
increases up to reaching a peak. As a rule, the impact of 
the increased motion cavity aspect ratio are smaller, but 
some exceptions occur. In fact, at any given aspect ratio 
larger than unity, should the Rayleigh number be such 
that heat transfer occurs by pure conduction, then speak- 
ing of motion intensity has no meaning and, obviously, the 
inclination of the cavity cannot have any effect on the heat 
transfer performance, which is for example what happens 
for A = 4 at Ra = 10°. On the other hand, should the Ray- 
leigh number be of the order of that corresponding to the 
onset of convection, then the inclination of the cavity 
may result in a degradation of the heat transfer perfor- 
mance, due to the fact that at large tilting angles conduc- 
tion tends to become the dominant heat transfer mode, 
which is for example what happens for A=4 at Ra=10”. 


Wmax = 26.97 Wl max = 24.94 


Fig.6 Isotherm and streamline contours plots for Ra = 10°, A= 4 and @ = 0° — 60° 
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Fig. 7 Distributions of Nu vs. Ra for ọ = 0° and A = 0.25 — 4 
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Fig.8 Distribution of Nu vs. A for ọ = 0° and Ra = 10° — 107 
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Fig. 9 Distributions of Nu/Nug vs. ọ for A = 0.25 and Ra = 
10*— 107 
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Fig. 10 Distributions of Nu/Nuy vs. @ for A= 1 and Ra = 10° 


= 10! 
1.1 
A=4 
= 
Z 1.0 0-0 9-9 
E 
Z o Ra= 108 ~ 
a Ra= 105 
a Ra= 10° 
e Ra= 107 4 
0.9 pe a” a’ a’ o” o 
0 15 30 45 60 


Tilting angle ọ [deg] 


Fig. 11 Distributions of Nu/Nup vs. ọ for A = 4 and Ra = 10° 
— 107 


The value of ọ corresponding to the peak of (Nu/Nupo) 
is defined as the optimal tilting angle for maximum heat 
transfer and denoted as op. The value of the maximum 
average Nusselt number occurring at such optimal incli- 
nation is denoted as Nuopt- As the tilting angle of the en- 
closure is further increased above Pops, the ratio (Nu/Nup) 
decreases, which is due to the excessive growth of the 
mentioned stagnation regions that form on both the heated 
and cooled walls. Obviously, when the increased stagna- 
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tion effect outweighs the increased velocity effect, the 
ratio (Nu/Nuo) becomes lower than unity, which means 
that the inclination with respect to gravity leads to a de- 
gradation of the heat transfer performance of the cavity. 

As far as the dependence of the optimal tilting angle 
on the Rayleigh number and the aspect ratio of the cavity, 
a set of distributions of opt versus Ra are depicted in Fig. 
12 for different values of A, showing that, as long as 
convection is the dominant heat transfer mode, opt in- 
creases as Ra is decreased, and A is increased. The cor- 
responding distributions of Nuop: versus Ra are plotted in 
Fig. 13. 

After having outlined the basic heat and fluid flow 
features, and discussed the effects of A, Ra, and ọ, on the 
overall thermal performance of the enclosure, a number 
of dimensionless correlations are now proposed for pre- 
dicting the optimal tilting angle Popi and the correspond- 
ing average Nusselt number Numax, as well as the ratio 
(Nu/Nup) in the investigated range of tilting angles. 

The numerical results obtained for the Nusselt Num- 
ber for vertical cavities, i.e. ọ = 0°, for the whole set of 
Rayleigh number, with 10* < Ra < 10’, and aspect ratios, 
with 0.25 < A < 4, as far as convection heat transfer is 
prevalent, i.e. Nu = 2xA, may be correlated through the 
following algebraic equation, derived by a multiple re- 
gression procedure: 


45 
30 
g 

15} 

0 

10° 104 10° 108 10’ 
Rayleigh number Ra 
Fig. 12 Distribution of Pop vs. Ra 

100 


104 10° 10° 10’ 
Rayleigh number Ra 


Fig. 13 Distribution of Nug,; vs. Ra 
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Nuo = (0.142 + 0.068 - In(A))- Ra?2+0. 1A”? (10) 


with a 5.3% standard deviation of error and 92% of 
data within a +5.3% range of error. The numerical results 
obtained for other tilting angles, with 0° < ọ < 45°, in the 
same ranges of Ra and Nu, may be correlated through the 
following algebraic equation, derived by a multiple re- 
gression procedure: 


Nu dy, dp, P— 55: Ra? I6. A 
—— = (=. 06- Rae A se (TI 
where 
a, =-0.04+0.31-A°% (12) 
104 
Dra =a (13) 
Cra = 0.005- Ra??? (14) 


with a 2.4% standard deviation of error and 95% of 
data within a +4.3% range of error. 

The numerical results obtained for the optimum posi- 
tion tilting angle for the whole set of Rayleigh number, 
with 10‘ < Ra < 107, as far as convection heat transfer is 
prevalent, i.e. Nu => 2xA, may be correlated through the 
following algebraic equations, derived by a multiple re- 
gression procedure: 

for low aspect ratio cavities, 1.e.0.25<A<2 

0.5 <5 -0.22 
5-A Li (15) 
In(Ra- A ~~) 

with a 0.026 rad (1.5°) standard deviation and range of 
error from —0.055 rad (-3.1°) to +0.040 rad (+2.3°); 

for high aspect ratio cavities, i.e. 3 < A < 8 
(4.2 203 A) R04 (16) 

with a 0.029 rad (1.7°) standard deviation and range of 
error from —0.06 rad (-3.5°) to +0.02 rad (+1.3°). 

Absolute errors has been considered as the angle val- 
ues may be in a range between 0° to 45° that means that 
relative error of the correlation would be very dependent 
on the actual value, thus being less meaningful. 


opt = 0.155- A +30-Ra?® + 


Popt = 9.28 + 


Conclusions 


Natural convection in air-filled rectangular cavities in- 
clined with respect to gravity, so that the heated wall is 
facing upwards, has been studied numerically under the 
assumption of two-dimensional laminar flow, for height- 
to-width aspect ratios of the enclosure from 0.25 to 8, 
Rayleigh numbers based on the length of the heated and 
cooled walls from 10° to 10’, and tilting angles of the 
enclosure from 0° to 75°. A computational code based on 
the SIMPLE-C algorithm has been used for the solution 
of the system of the mass, momentum and energy trans- 
fer governing equations, assuming the initial condition of 
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motionless fluid at the uniform temperature of the cooled 
wall up to the achievement of a steady-state solution. 

It has been found that, as the inclination angle is in- 
creased starting from the vertical position, the heat trans- 
fer performance increases up to reaching a peak at an 
optimal tilting angle, that, in line of principle, unless the 
value of the Rayleigh number was that corresponding to 
the onset of convection or just higher, turned out to in- 
crease as the Rayleigh number was decreased, as well as 
the height-to-width aspect ratio of the cavity was in- 
creased. 

Dimensionless correlating equations have been pro- 
posed to evaluate the optimal tilting angle and the heat 
transfer performance of the enclosure. 
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